A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

03 (08.19.04; 08.19.5 SN 1994D; 11.09.1 NGC 4526; 10.07.2; 12.04.3) 



ASTRONOMY 

AND 
ASTROPHYSICS 



SN 1994D in NGC 4526: a normally bright type la supernova 



on 
on 

ON 



Georg Drenkhahn 1 2 and Tom Richtler 1 

1 Sternwarte der Universitat Bonn, Auf dem Hiigel 71, 53121 Bonn, Germany (richtler@astro.uni-bonn.de) 

2 Max-Planck-Institut fur Astrophysik, Postfach 1523, 85740 Garching bei Miinchen, Germany (georg@mpa-garching.mpg.de) 

Received 17 March 1999 / Accepted 19 July 1999 



Of 

m : 
t> : 

CO ■ 

> ■ 
in : 

CO ■ 

r» ■ 

o : 

ON . 

ON ■ 



Of 
i 1 

o 



X 



Abstract. SN 1994D of type la has been suspected not to fit 
into the relation between decline rate, colour, and brightness. 
However, an individual distance of its host galaxy, NGC 4526, 
other than that of the Virgo cluster, has not yet been pub- 
lished. We determined the distance by the method of globular 
cluster luminosity functions on the basis of HST archive data. 
A maximum-likelihood fit returns apparent turn-over magni- 
tudes of 23.16 ± 0.16mag in V and 21.96 ± 0.09mag in 
I. The corresponding distance modulus is 30.4 ± 0.3 mag, 
where the error reflects our estimation of the absolute distance 
scale. The absolute magnitudes (not corrected for decline rate 
and colour) are —18.67 ± 0.30 mag, —18.62 ± 0.30 mag, and 
— 18.40±0.30mag for B, V, and I, respectively. The corrected 
magnitudes are —18.69 ± 0.31 mag, —18.69 ± 0.31 mag, and 
-18.44 ± 0.31 mag. 

Compared with other supernovae with reliably determined 
distances, SN 1994D fits within the errors. It is therefore 
not a counter-example against a uniform decline-rate-colour- 
brightness relation. 

Key words: supernovae: general - supernovae: individual: 
SN 1994D - galaxies: individual: NGC 4526 - globular clus- 
ters: general - distance scale 



(1995), based on a comparison with SN 1989B in M66 and 



1. Introduction 

Supernovae (SNe) of type la are the best standard candles in 
the universe. Hamuy et al. ( 1996 » showed how remarkably 
small the brightness dispersion among SNe la is, for a sam- 
ple of high observational quality. Moreover, decline rate and 
colour at maximum phase can be used to introduce correc- 
tions and to squeeze the dispersion even below 0.15 mag in 



B, V, and I (Tripp 1998; Richtler & Drenkhahn 1999). But 



a few SNe have been observed which apparently do not fit 
with the Calan/Tololo sample of Hamuy et al.; the prominent 
cases being SN 1991bg in NGC 4347, SN 199 IT in NGC 4527, 
and SN 1994D in NGC 4526, which are too bright/faint/red etc. 
However, checking the literature, one finds that the claim for 
anomaly often relies on indirect arguments, i.e. without a reli- 
able knowledge of the distance of the host galaxy. For exam- 
ple, SN 1994D was labelled "anomalous" by Richmond et al. 



SN 1980N in NGC 1316, which made SN 1994D appear too 
bright by 0.5 mag. A disadvantages is that M66 only has a 
Tully-Fisher distance, for which it is known that large resid- 
uals to the mean relation exist for individual galaxies. More- 
over, this method suffers from uncertain extinction, which, to- 
gether with the dispersion in the Tully-Fisher relation, results 
in an error in the ab solute magnitude of the order 0.25 mag 
(Jacoby et al. 1992 ). On the other hand, the distance mod- 
ulus of NGC 1316 was adopted as 31.02 mag while globular 
clusters (GCs), Cepheids, and s urface brightness fluct uation 
(SBF) measu remen ts (Kohle et a l.|l996| Delia Valle et al. |l998 , 
Madore et al. 1998, Jensen et al. 1998) now agree on a value of 



~ 31.35 mag for the Fornax Cluster distance modulus. Due to 
the cluster's compactness we identify the distance to NGC 1316 
with the cluster distance. An up-to-date discussion of the For- 
nax Cluster distance can be found in Richtler et al. ( 1999 ). 

Sandage & Tammann ( |1995[ ) also found SN 1994D too 
bright by 0.25 mag. They used the fact that NGC 4526 is a 
member of the Virgo cluster and based their comparison on 
other Virgo SNe. While the difference may be marginal, any 
correction for the decline rate would increase the discrepancy, 
hence rendering the decline-rate-luminosity relation not gen- 
erally valid. However, no individual distance for NGC 4526 
has yet been published. The distance value used by Richmond 
(1995) and Patat (1996) is an unpublished SBF distance modu- 
lus of 30.68 ±0.13 mag from Tonry. In this paper, we shall give 
a distance for NGC 4526, derived by the method of g lobula r 
cluster luminosity functions (GCLFs) (e.g. Whitmore 1997 ), 
and we want to demonstrate that the brightness of SN 1994D 
fits reasonably well to the Fornax SNe. 



2. Data and reduction 

2.1. Data 

Our data consist of HST-WFPC2 images and have been taken 
from the HST archive of the European Coordinating facility at 
ESO. Table [j] lists the exposure times and filters associated with 
the data. The date of the observation was May 8th, 1994 (Ru- 
bin 1994[ ). Pipe-line calibration was done with the most r ecent 
calibration parameters according to Biretta et al. ( 1996 ). All 
frames were taken with the same pointing and the photometry 
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Table 1. The names of the respective HST data set, the filter, 
and the exposure time in seconds. 



Name 


Filter 


^exp 


u2dt0501t 


F555W 


60 


u2dt0502t 


F555W 


230 


u2dt0503t 


F555W 


230 


u2dt0504t 


F814W 


60 


u2dt0505t 


F814W 


230 


u2dt0506t 


F814W 


230 



was performed on averaged frames with an effective exposure 
time of 520 s each. 

Hot pixels were removed with the IRAF task WARMPIX. 
The limited Charge Transfer Efficiency was corrected for with 
the prescriptions given by Whitmore & Heyer (1997 ). Cosmics 
were identified with the IRAF task CRREJ and eventually the 
frames were combined with a weighted average. 

We modelled the galaxy light with the IRAF task IMSUR- 
FIT and achieved a very satisfying subtraction of the galaxy 
light with no visible residuals. After this, the photometric ef- 
fects of the image distortion were corrected fo r by a multipli- 
cation with a correcting image (Holtzman et al. 1995). 



2.2. Object search and photometry 

Due to the Poisson noise of the galaxy light, one finds a 
noise gradient towards the centre of the galaxy. The normally 
used finding routin e DAOFIND from the DAOPHOT pack- 
age (Stetson 1987) does not account for a position-dependent 
background noise. Using DAOFIND with a mean noise level 
results in many spurious detections near the galaxy centre and 
missed objects in regions with lower noise background. There- 
fore we used the fi nding routine of the SEXTRACTOR program 
(Bertin & Arnouts 1996 ). The lower threshold was adjusted in 
order to include also the strongest noise peaks with a number 
of connected pixels of 3 and a threshold of 1.4 • a. Our ini- 
tial object list consists of 304 found and matched objects in the 
three wide-field frames of the F555W and F814W images. The 
planetary camera data were not used because they showed only 
very few globular cluster-like objects. 

For these objects aperture photometry (with DAOPHOT) 
has been obtained through an aperture of 0'.'5. These instrumen- 
tal magnitudes were finally corrected for total flux and trans- 
formed into standard Johnson magnitudes with the formulae 
given by Holtzman et al. (1995). 



3. The luminosity function 

3.1. The selection of globular cluster candidates by 
elongation, size and colour 

The selection of cluster candidates in HST images benefits 
greatly from the enhanced resolution. While in ground-based 
imaging, cluster candidates often cannot be distinguished from 
background galaxies, this problem is less severe in HST im- 
ages. Galaxies that are so distant that their angular diameter is 



comparable to cluster candidates in NGC 4526, normally have 
quite red colours. Remaining candidates are foreground stars 
and faint background galaxies, where the redshift balances an 
intrinsic blue colour. Our strategy, described in more detail be- 
low, is the following: we first select our sample with respect to 
the elongation as given by the S EXTRACTOR program, then the 
angular size measured as the difference between two apertures, 
and finally the colour. 

3.1.1. Elongation 

The elongation is taken as the ratio E = A/B, where A and B 
are the major and the minor axes as given by SEXTRACTOR. 
However, the undersampling can cause a spurious elongation if 
a given object is centred between 2 pixels. We therefore chose 
E < 1.9 in both filters to exclude only the definitely elongated 
objects. 

3.1.2. Size 

The SEXTRACTOR stellarity index is not very suitable for 
undersampled images because the neural network of SEX- 
TRACTOR is trained with normally sampled objects (Bertin & 
Arnouts 1996). We use the magnitude difference in two aper- 
tures as a size measure (Holtzman et al. 1996| ). These aper- 
tures should be compatible with the expected size of a GC in 
NGC 4526. A typical half-light radius of a GC is of the order of 
a few pc. Thus the larger clusters are marginally resolvable, as- 
suming a distance of 13 Mpc to NGC 4526, which corresponds 
to 6 pc per pixel. The parameter mi2 := m,\ — m^— (mi- 7712)* 
is defined as the magnitude difference between one and two 
pixel apertures mi — m,2 minus the appropriate value for a star- 
like object (7711 — 7712)*. The latter subtraction accounts for the 
position-dependent PSF of the WFPC2. 

To find a reasonable size limit, we simulated cluster can- 
didates by projecting the data for Galactic GCs at a distance 
of 13 Mpc, based on the McMaster catalogue of Galactic glob- 
ular clusters (Harris 1996). We generated oversampled King- 
profiles (King 1962) with the geometric parameters from this 
catalogue and convolved them with an oversampled PSF of the 
WFPC2. The resulting images were rebinned at four differ- 
ent subpixel positions giving WFPC2 images of the Galactic 
globular cluster system (GCS) as it would be observed from a 
distance of 13 Mpc. Fig. [l] shows the measured m\i values as 
functions of the core radii r c . It turns out that TO12 is a suffi- 
ciently good quantity to measure r c . Applying a selection cri- 
terion of TO12 < 0.45 in both filters excludes extended GCs 
with r c > 5pc. This reduces the size of the GC sample but 
does not introduce a bias because r c and My are not corre- 
lated. The mi2 magnitude is not suitable to place a lower size 
limit because the smaller GCs are unresolved and cannot be 
distinguished from a stellar-like image. 

Fig. |2| compares the distribution of 777,12 of the sample of 
our detected sources with that of the simulated clusters for both 
filters. The dotted histogram represents the background sources 
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Fig. 1. The measured ran values of the simulated Galactic GCs 
as a function of the core radius r c . Every GC is represented 
by 12 crosses according to the 3 WF cameras and 4 different 
subpixel positions. 



(see Sec. 3.3). The objects with to 12 > 0.45 can be almost 
entirely explained as resolved background objects. 

3.1.3. Colours 

The range of existing colours among Galactic GCs is not a 
very good guide in our case: large error bars and a possibly 
different metallicity distribution may cause a different shape 
of the colour histogram. However, as Fig. |3| shows, a selection 
with regard to shape and size selects also according to colours, 
i.e. the redder objects, which are presumably redshifted back- 
ground galaxies, can be effectively removed. 

The colour histogram of the remaining objects shows some 
similarity to other GCSs with regard to the apparent "bimodal- 
ity", for example in M87 (Whitmore et al. 1995 ). One may sus- 
pect that there is a dichotomy in the population of GCs in the 
sense that the blue (=metal poor) objects belong to a halo pop- 
ulation and the red ones to the bulge. Unfortunately, the num- 
ber of objects found does not allow a deeper investigation. The 
mean colour is (V — I) = 1.1. 

We also exclude very faint objects with large photomet- 
ric errors Am > 0.3 in at least one of the two filters. They 
hardly contribute to the determination of the turn-over magni- 
tude (TOM) because of their large uncertainty. 



3.2. Completeness 

The completeness /(to) ■ dm is the probability of finding an 
object in the magnitude interval [to, TO+d77i] . On frames with a 
homogeneous background, it depends mainly on the brightness 
of an object. In our case, where the object search was done on 
a frame with a subtracted galaxy, the noise of the background 
p varies strongly, and must be treated as a second parameter in 
I(m;p). 
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Fig. 2. These histograms show the distribution of TO12 for 
the detected sources (upper panel) and the simulated clusters 
(lower panel). The dotted histogram shows the expected back- 
ground objects (=galaxies) derived from the Medium Deep Sur- 
vey frames. Objects with values TO12 > 0.45 are resolved and 
mostly background objects. The distribution for F814W is es- 
sentially the same. 



We evaluated the completeness in our V, I sample by in- 
serting artificial objects in the V-frame and we used the same 
coordinates for inserting objects in the /-frame, after having 
fixed the colour for all objects to the value of V — I = 1.1. 
In total we used 19000 artificial "stars" in the magnitude range 
23 < V < 25.5 and evaluated the completeness in bins of 
0.1 mag. The noise as the second parameter is taken to be the 
standard deviation p of pixel values in the annulus used by the 
aperture photometry. In practice, we used three completeness 
functions for three different p-values 1.75 DN, 2.05 DN, and 
3.35 DN, which represent the intervals [1.6, 1.9], [1.9, 2.2], and 
[2.2, 4.5], respectively. 

Fig. Q shows the completeness for the three noise intervals 
and for the two filters. The distinct difference between the three 
curves demonstrates that p is indeed a reasonable parametrisa- 
tion. 



4 



Georg Drenkhahn & Tom Richtler: SN 1994D in NGC 4526: a normally bright type la supernova 



35 - 

30 

25 - 

I 20 
E 

* 15 f- 
10 

5 - 




-0.5 




1.5 2 2.5 



Fig. 3. This colour histogram shows the shape- and size- 
selected sample of objects. The dotted line denotes the distribu- 
tion of the complete sample. It can be seen that the shape and 
size distribution can already discriminate against red objects 
which are presumably distant galaxies. For the final sample, 
we chose the colour interval 0.8 < V — I < 1.4. 
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Fig. 4. The completeness for V and / for three different noise 
intervals, which are characterised by the mean values p = 
1.75 DN, 2.05 DN and 3.35 DN (from right to left) in the V 
filter. 



Table 2. The exposures from the Medium Sky Survey which 
have been used for the background evaluation. Given are the 
data set names, total exposure times in seconds, and the number 
of individual frames, which entered the averaging process. 





^cxp [s] 


^cxp [s] 


# 


# 


Name 


in F606W 


in F814W 


in F606W 


in F814W 


u30h3 


1200 


1266 


6 


8 


uy400 


866 


1000 


6 


6 


uy401 


600 


700 


4 


8 


uyxlO 


800 


700 


3 


3 


uzxOl 


696 


1550 


5 


4 



3.3. Foreground and background sources 

The field of our HST images is too small for a proper analy- 
sis of the contaminating background sources that survived all 
selection processes. 

Therefore, we used t he W FPC2 Medium Deep Survey 
(MDS) (Ratnatunga et al. 1999 ) to obtain an estimation of the 
background. Suitable fields should be possibly nearby, in par- 
ticular regarding the Galactic latitude, and should be deep with 
several single exposures. We found the exposures which are 
listed in Table || together with some information. The frames 
are deep enough (in F606W and F8 14W) that the source finding 
is complete in the relevant magnitude interval. V and I magni- 
tudes are calculated according to the prescriptions of Holtzman 
et al. (119951). 



The selection criteria discussed above are applied in exactly 
the same way as in the case of the NGC 4526 frames. It turns 
out that most background objects are selected out. The Poisson 
scatter of background objects from frame to frame is in the ex- 
pected statistical range. Thus any clustering in the background 
galaxy number density does not influence the background LF. 
So we can assume that the background LF derived from the 
MDS frames gives also a good estimate for the background LF 
near NGC 4526. 

As can be assessed also from Fig. [z[ the fraction of back- 
ground objects in the selected samples is only 7%. The result- 
ing luminosity function b(m) enters the maximum-likelihood 



fitting process as described in Sec. 3.4.1 



3.4. The turn-over magnitude 

3.4.1. Remarks on the maximum-likelihood fit and the TOM 
of NGC 4526 

For the representation of the GCLF, we choose the t$ -function, 
which has been shown to be the best analytica l repr esentation 
of the Galactic GCLF and that of M31 (Seeker \l99%. It reads 



t 5 (m; m°,(T t ) = 



3V5- 



(m — m ) 2 
~5aJ 



(1) 
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The maximum-likelihood fitting method evaluates the most 
probable parameters for a given distribution <f> by maximising 
the likelihood 



L(m°,(Jt) = W_4>{m l \m a ,a t ). 



(2) 



For the the oretic al backgroun d of th is metho d see B evington & 
Robinson (]1992| ), Caso et al. ( |l998| ). Seeker ( |l992| ) and Seeker 
& Harris ( |1993| ) already applied it to the fitting of GCLFs. 

What is to be fitted is the t§ -function modified by the differ- 
ent influencing factors, i.e. the completeness, the background, 
and the photometric error s. We used an extended version of 
the distribution of Seeker ( 1992| ) that accounts for a position- 
dependent background noise. 



<p(m;p,m°,a t ) = 

K ■ k(to; p, to , a t ) ± (1 — K) ■ f3{m\ p) 

with k being the normalised distribution of the GCs 

n(m; p,m° ,a t ) = 

/(to; p) ■ (/ tc,(m'; to , a t ) ■ s(m; to', p) dm.') 



(3) 



(4) 



and (3 being the normalised distribution of the background ob- 
jects 



/3(to; p) = /(to; p) • b(m). 



(5) 



Special care must be taken when n and j3 are normalised and 
mixed through the mixing value K, because this depends on 
the completeness and therefore on p. 

e is a Gaussian distribution accounting for the photometric 
error 



e(to; to', p) 



exp 



(to — m') 



IY2 



/2tt • ct(to',p) * \2-a(m',p) 
with a magnitude- and noise-dependent width 

cr(m, p) = y/lQOMm-a) + p 2 . 10 0.8-(m-6) 



(6) 



(7) 



as used by DAOPHOT (Stetson |1987| ). a and b depend on the 
photometry parameters but are more easily determined by just 
fitting the photometry data to (Q). 

The necessity for this more complicated approach is most 
clearly seen in the noise-dependent completeness functions in 
Fig. fj| One could also perform the analysis for objects in a 
small range of background noise but this would further reduce 
our already small GC sample. 

Fig. |] shows two-dimensional slices through the three- 
dimensionally fitted luminosity function at a noise level of 
p = 1.75, 2.00, and 3.25 DN, respectively, in both filters. The 
background noise of p = 2.00 DN corresponds approximately 
to the mode of the background noise distribution and is there- 
fore most representative. The figure also displays the low back- 
ground contribution and a down-scaled completeness function 
for p = 2.00 DN. The histogram includes all objects indepen- 
dent of their p- values. It therefore does not correspond to a sin- 
gle fitted function but is rather a superposition of many func- 
tions belonging to a variety of p-values. Though the histogram 



Table 3. The fitted turn-over magnitudes and tg dispersions for 
the globular cluster system of NGC 4526. 



Filter 



mo 



V 

I 



23.16 ±0.16 
21.96 ±0.09 



1.01 ±0.12 
0.94 ± 0.09 



itself is not used in the fitting it shows that the fitted distribution 
models the data reasonably well. From simply looking at the 
histogram one could conclude that the TOM lies at a slightly 
fainter magnitude. But one has to keep in mind the large statis- 
tical scatter within the single bins and that this picture changes 
with different bin positions and widths. Note also the differ- 
ences between the maxima of the functions 4>(m; 1.75 DN), 
0(to;2.OODN), and 0(to; 3.25DN) in V caused by the dif- 
ferent completeness limits. The completeness function in our 
data decreases at the same magnitude value as the declining 
part of the intrinsic LF. It is therefore crucial to determine the 
completeness function most precisely and any error in this part 
of the calculation introduces a significant systematic error in 
the TOM. This problem is only resolved if the observation is 
complete down to at least one magnitude beyond the TOM. 

Fig. ^ shows the likelihood contours for the fitted param- 
eters and Table [| gives the resulting values. The width at is 
0.09 mag larger than at for the Galactic system (Table ^) but 
still small compared to other early-type galaxies (Whitmore 



1997 use at = 0.78 • a g to convert the Gaussian widths). 



3.4.2. The TOM of the Galactic system 

Although the Galactic TOM has been discussed by many au- 
thors, it is not a parameter that can be simply looked up in the 
literature. The exact value depends on the method of fitting, on 
the selection of GCs that enter the fitting process, on the dis- 
tance scale of GCs, and on the exact application of a reddening 
law. Even if our values do not differ much from previous eval- 
uations, we discuss them for numerical consistency. 

We take the LMC distance as the fundamental distance for 
calibrating the zero-point in the relation between metallicity 
and horizontal branch/RRLyrae brightness. The third funda- 
mental distance determination beside trigonometric parallaxes 
and stellar stream parallaxes is the method of Baade-Wesselink 
parallaxes. It has been applied to the LMC in its modified form 
known as the Barnes-Evans method. So far, it has been applied 
to Cepheids in NGC 1866 (Gieren et al. |1994| ), and the most 
accurate LMC distance determination to date stems from the 
period-luminosity relation of LMC Cepheids by Gieren et al. 
(1998). We adopt the distance modulus from the latter work, 
which is 18.46 ± 0.06 mag, and which is in very good agree- 



ment with most other work (e.g. Tanvir 1996 ) 



If we adopt the apparent magnitude of RR Lyrae stars in the 
LMC from Walker fll992[ ), 18.94 ± 0.1 mag for a metallicity 
of [Fc/H] = — 1.9dex, and the metallicity dependence from 
Carretta et al. ( |1999[ ), we find 



M V {RR) = (0.18 ± 0.09)([Fe/H] ± 1.6) ± 0.53 ± 0.12. (8) 
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Fig. 5. The result of the maximum-likelihood fit for the 
V and / filter (solid line). Histogram: all selected objects 
independent of their p values. Solid lines: fitted GCLFs 
</>(m;1.75DN), 0(m;2.OODN), and </>(m; 3.25DN). Dotted- 
dashed line: down-scaled completeness 0.6 • J(m; 2 DN). 
Dashed line: background contribution (3(m; 2DN). The ordi- 
nate is the dimensionless probability per magnitude for (f>, j3 
and the histogram while for / it just denotes the probability. As 
required by any probability density function the area under the 
fitted distribution function (f>(m; p) and the histogram is unity 
for all values of p. 



This zero-point is in excellent agreement with the value of 
0.58 ± 0.12 mag derived from HB-brightnesses of old LMC 



globular clusters from Suntzeff et al. (1992), if the above metal 



licity dependence is used. It also agrees within a one-cr limit 



with 0.56 ± 0. 12 ma g (Chaboyer |1999|) and 0.57 ± 0.04 mag 
(Carretta et al. 



1999) 



We calculate the absolute magnitudes of Galactic GCs from 



the McMaster catalogu e (Har ris |1996[ ), using Ay /E(B-V) = 

-I)/E{B-V) 



3.1(Rieke&Lebofsky 
(Ardeberg & Virdefors 



1985) and EjV- 



1.53 



1982) as a compromise between the ex- 
treme values 1.6 mag (Rieke & Lebofsky |i985] ) and 1.35 mag 



(Drukier et al. |1993j). We exclude clusters with E(B - V) > 1 
so that the sample embraces 120 clusters in V and 93 in I. 
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Fig. 6. 1, 2, and 3er contours for the fitted parameters m° and 
at in V and /. 

Table 4. Fitted TOMs and a t -values for the Galactic GC sys- 
tem. The / values are biased due to selection effects as dis- 
cussed in the text. 



Filter 



M° 



at 



# of GCs 



V 
I 



-7.61 ±0.08 
-8.60 ±0.07 



0.92 ±0.07 
0.85 ± 0.06 



120 

93 



Fig. ^Jshows the scaled histograms (V and /) for the sample 
of Galactic globular clusters. The solid line is the t$ -function 
resulting from the maximum-likelihood fit. It is not the fit to the 
histograms, which are shown just for visualisation purposes. 

It is apparent from Fig. ^that the symmetry in the V distri- 
bution breaks down for clusters at the faint end of the distribu- 
tion. It is reasonable to include in the fit only the symmetric part 
up to a certain limiting magnitude. To determine an appropriate 
cut-off magnitude m c we carry out the maximum-likelihood fit 
for various m c and look at how the fitted parameters evolve. 
Fig. H shows that the V-TOM M v and the width a t clearly de- 
pend on m c . The inclusion of the asymmetric part shifts My 
systematically to greater values. This trend stops at our adopted 
value m c = — 5.8 mag and lowering m c further just results in a 
random scatter of My. 

There is no asymmetric part in the histogram of the / band 
data in Fig. The reason for this is a selection effect in the 
/ photometry. Looking at the McMaster catalogue data reveals 
missing / data in approximately 22% of the listed GCs. All 
of these clusters with missing data belong to the fainter end 
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Fig. 7. This graph shows the normalised distribution of ab- 
solute magnitudes in V and / for Galactic globular clusters. 
Overplotted as solid lines are the is-functions resulting from a 
maximum-likelihood fit. The bold histogram in the left panel 
shows the GCs with missing /-band photometry. The dotted 
vertical line denotes the cut-off magnitude at rn c = —5.8. 
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Fig. 8. Dependence of the fitted M v (symbol X) and at (sym- 
bol +) on the cut-off magnitude m c for the Galactic GCS. The 
chosen m c — —5.8 is indicated by the dashed line. 



of the luminosity function in V. This is shown by plotting the 
subset of GCs without / magnitudes onto the My histogram in 
Fig. ^. This incompleteness shifts the fitted / TOM to brighter 
magnitudes. Thus, a comparison between the difference of the 
fitted TOMs and the mean colour of the Galactic GCs gives 

0.87 ± 0.09 = (V - 7) MW < M v - M° = 1.01 ± 0.11. (9) 

The median is equal to the mean of the V — I distribution for 
the Galactic GCs. We therefore use 



My 



(V-I) 



MW 



.48 ±0.10 



(10) 



as an absolute TOM in the / band. Alternatively, we intro- 
duce a cut-off magnitude in the / band data and fit only the 
more complete luminous part of the distribution like we did in 
the V band. As expected, the fitted / TOM becomes smaller 
with decreasing cut-off magnitude. If the cut-off magnitude 
approaches the TOM the fitted TOM values begin to scatter 
greatly around approximately —8.52 mag because the sample 
size decreases and the fit results are less well determined. We 
therefore do not rely on this approach but rather regard it as a 
confirmation of the validity of (^0[). 

3.5. Metallicity correction and distance 

The luminosity of a GC depends not only on its mass but also 
on its metallicity. Any difference between the metallicities of 
the GCS of NGC 4526 and the Galactic GCS must be corrected 
to avoid a systematic bias. 

Th e amo unt of this correction has been studied by Ashman 
et al. (1995 1 and we adopt their values. The only metallicity 
indicator that is available for us is the mean V — I colour of the 
GCS, for which we adopt the relation given by Couture et al. 
<\l99ty: 



V - I = 0.2 • [Fc/H] + 1.2. 



(11) 



A mean colour of (V — I) = 1.1 indicates [Fe/H]NGC4526 = 
—0.5 ± 0.3 and we find a metallicity correction of AM V = 
0.28 ± 0.19 and AM? = 0.11 ± 0.05. The large error of 
±0.3 dex for the metallicity arises from the fact that (V — I) 
depends only weakly on the metallicity. This value can be esti- 
mated from looking at the scatter in the [Pe/H]— (V — /) plot 
of the Galactic GCS. 

Hence the distance moduli in V and / are 



Hv = 30.49 ± 0.26 and m = 30.33 ± 0.14. 



(12) 



If the universality of the TOM of ps 0.2 mag (Whitmore [1997p 
and the uncertainty in the horizontal branch calibration of « 
0.1 mag is included as external errors we obtain 

AtNGC4526 = 30.4 ± 0.3 mag (13) 

and 

£>NGC4526 = 12.0±1.6M P C (14) 

as our final value for the distance modulus and distance. 
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Table 5. The table lists the coefficients according to ( |l5[ ) for 
the relation between brightness, decline rate, and colour for the 
Calan/Tololo sample. 



-1.5 



Band 



R 



B 
V 
I 



-3.306 ± 0.063 
-3.320 ±0.058 
-3.041 ± 0.060 



-0.48 ±0.23 
-0.52 ±0.19 
-0.40 ±0.22 



-1.51 ±0.62 
-0.83 ±0.57 
-0.81 ±0.59 



Thus, NGC 4526 is located in the foreground of the Virgo 
cluster and a comparison with apparent magnitudes of other 
SNe in the Virgo cluster is misleading. This result is consis- 
tent within la- with Tonry's unpublished SBF distance modulus 
30.68 ± 0.13mag, which was used by Richmond et al. (1995). 



4. Comparison with other supernovae 

4.1. The dependence on decline rate and colour 

The publication of the Calan/Tololo-sample of la SNe (Hamuy 
et al. 1996) confirmed in a compelling way what has been sug- 
gested earlier, namely that the luminosity of la's depends on 
the decline rate and also on the colour of the SN at maximum 
phase. 

Here we report convenient expressions which relate in the 
sample of Hamuy et al. the decline rate Amis and the fore- 
ground extinction-corrected colour i? max — V maK with the ap- 
parent maximum magnitude (the recently published sample by 
Riess et al. (1998) does not change significantly the numerical 
values). These linear relations read 



+ b ■ (Am 15 
= 5 • log cz ± Z 



1.1) + R ■ (B n 



(15) 



where i E {B, V, 1} denotes the different bands. Table |5] lists 
the coefficients, which have been obtained by a maximum- 
likelihood fit, together with their errors. Only SNe with i? max — 
Knax < 0.2 have been used for the fit. Fig. |] shows the relation 
between B max — VJ nax and the apparent maximum brightness 
corrected for decline rate. The fact that the red SNe are also 
well fitted by this relation indicates the validity of this relation 
for a larger colour range. The coefficients of the colour term 
also deviate strongly from the Galactic reddening law so that 
one may conclude that the red colours are indeed intrinsic and 
not significantly affected by reddening. 

The corrected magnitudes can now be compared with those 
of other SNe. 

4.2. Comparison with other SNe 

A meaningful comparison of SN 1994D with other SNe re- 
quires the distances of the comparison objects to be derived 
with the same method of GCLFs, because then the comparison 
is valid without the need for an absolute calibration of GCLFs. 
Moreover, these SNe should be well observed. These two con- 
ditions are fulfilled only in the Fornax galaxy cluster with 



E 

< 

+ 




Fig. 9. This plot shows the relation between B max — V max rate 
and the apparent maximum brightness in B corrected for the 
decline rate. A clear dependence is visible. We point out that 
only blue SNe with i? m ax — Vmax < 0.2 have been included in 
the fit. The fact that the redder SNe are well represented by this 
relation is a point in favour of the validity of this relation even 
over a larger colour range. 



its 3 la SNe SN 1980N, SN 198 ID (both in NGC 1316), and 
SN 1992A (in NGC 1380). The light c urves o f the first two SNe 
have been presented by Hamuy et al. (1991 ), while SN 1992 A 
(one of the best observed SNe) is described by Suntzeff ( 1996 ). 

All these SNe are relatively fast decliners, and thus fol- 
low the trend already visible in the Calan/Tololo SN sample 
(Hamuy et al. 1996| ) that spiral galaxies host all kinds of SNe, 
while early-type galaxies tend to host only fast decliners (which 
are intrinsically fainter). 

The Fornax cluster is well suited to serve as a distance stan- 
dard. It is compact with hardly any substructure and consists in 
its core region almost entirely of ellipticals and SO galaxies, 
the characteristics of a well relaxed galaxy cluster. The dis- 
tance based on GCLFs has been determined through the work 
of Kohle et al. ( |1996| ) and Kissler-Patig et al. ( |l997a| ), who 
analysed the GCSs for several early-type galaxies. In addition, 
Kissler-Patig et al. ( |l997b| ) and Delia Valle et al. ( |1998| ) deter- 
mined the luminosity function for the GCS of NGC 1380, host 
of SN 1992 A. Kohle et al. ( |1996[ ) derived mean apparent TOMs 
for their sample of 5 Fornax galaxies of V° = 23.67 ± 0.06. 
This is in perfect agreement with the TOM of NGC 1380 given 
by Delia Valle et al. ([1998]), which is V° = 23.67 ±0.11. The 



difference in distance moduli between Kohle et al. (1996) and 
Delia Valle et al. (1998) stems from a revised absolute TOM 



of the Galactic GCS. If this is based on the Gratton et al. 
(1997) relation, which is similar to our eq. (Q), the distance 
modulus is (iF om ax = 31.35 ± 0.15 mag. It is very satisfac- 
tory that several other independent methods confirm our dis- 



tance estimate (Richtler et al. 1999). For example, Madore et al 



( |1998| ) recently published a Cepheid distance to NGC 1365 
and quoted 31.35 ± 0.2 mag. Moreover, the method of sur- 
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Table 6. The observational data and distance moduli of the SNe. The first three rows show the extinction-corrected apparent 
magnitudes computed with A B /E{B - V) = 4.09, A V /E(B - V) = 3.09 and Aj/E{B - V) = 1.49 (Rieke & Lebofsky 



1985). Photometric data, colour excesses and decline rates are taken from Richmond et al. (1995) for SN 1994D, from Hamuy 



et al. (1996) for SN 1992A and SN 1980N, and from Hamuy et al. (1991 ) for SN 1981D. The absolute magnitudes are corrected 



for decline rate and colour by (1151) 





SN 1994D 


SN 1992A 


SN 1980N 


SN1981D 


^raax 


11.73 ± 0.02 


12.57 ±0.03 


12.49 ±0.03 


12.59 


Vmax 


11.78 ± 0.02 


12.55 ±0.03 


12.44 ± 0.03 


12.40 


^max 


12.00 ±0.05 


12.80 ±0.03 


12.70 ±0.04 




Am 15 (B) 


1.31 ±0.08 


1.47 ±0.05 


1.28 ±0.04 




^max Vmax 


-0.05 ± 0.04 


0.02 ± 0.04 


0.05 ± 0.04 


0.19 


E(B - V) 


0.04 ±0.03 








0(7) 




30.40 ± 0.30 


31.35 ±0.15 


31.35 ±0.15 


31.35 ±0.15 


M B , cor 


-18.69 ±0.31 


-18.99 ±0.16 


-19.02 ±0.16 




My :COT 


-18.69 ±0.31 


-19.01 ±0.16 


-19.04 ±0.16 




Mr, cor 


-18.44 ±0.31 


-18.69 ±0.16 


-18.74 ±0.16 





face brightness fluctuations gives a For nax distance modulus 
of 31.32 ± 0.24 mag (Jensen et al. |l998| ). 

For the other host galaxy in Fornax, N GC 13 16, no pub- 
lished GCLF distance exists. Grillmair et al. (1999), using HST 
data, found no turn-over but an exponential shape of the GCLF. 
They explained this with a population of open clusters mixed 
in, presumably dating from the past merger event. Preliminary 
results from ongoing work (Gomez et al., in preparation), how- 
ever, place the TOM of the entire system at the expected mag- 
nitude, but indicate a radial dependence in the sense that the 
TOM becomes fainter at larger radii. At this moment, we have 
no good explanation. Perhaps, the different findings correspond 
to the different regions in NGC 1316 investigated. This issue is 
not yet settled. 

An individual distance value for NGC 1316 based on plan- 
etary nebulae seems to indicate a shorter distan ce to this 
galaxy than to the Fornax cluster. McMillan et al. ( 1993 ) find 
Mpn, ngc 1316 = 31.13mag, but this value is based on a 
M3 1 distance modulus of 24.26 mag. Adopting the now widely 
accepted M31 Cepheid distance modulus of 24.47 mag (e.g. 
Stanek & Garnavich |1998|, Kochanek |1997|), the NGC 1316 



modulus increases to 31.34 mag, perfectly consistent with the 
global Fornax distance. We therefore adopt the Fornax distance 
modulus of /iFornax = 31.35 ± 0.15mag for SN 1992A and 
SN 1980N. 

Table ^ lists the absolute magnitudes of SN 1994D together 
with two Fornax cluster SNe (SN 198 ID was badly observed), 
corrected for decline-rate and colour. All values can be derived 
from the photometric data and distances shown in table ^ to- 
gether with our correction equation (|i"5|). 

One notes that SN 1994D seems to be 0.3 mag too faint in 
all three photometric bands. However, the error bars overlap 
with those of the other SNe. Thus, there is no reason to claim 
that this SN is peculiarly dim. On the other hand, if one takes 
Tonry's unpublished SBF distance modulus, which is 0.28 mag 
greater than ours, all three SNe brightnesses would agree 
within ±0.05 mag. At any rate, it is apparent that SN 1994D is 
definitely not overluminous and cannot be taken for a counter- 



example against the general validity of the decline-rate-colour- 
luminosity relation. 

4.3. The Hubble constant 

Once the absolute brightness of a supe rnova is known, the Hub- 
ble diagram of SNe la (Hamuy et al. 1996) allows the deriva- 
tion of the Hubble constant with an accuracy which is largely 
determined only by the error in the supernova brightness. The 
fact that the Fornax SNe are well observed in combination with 
a reliable distance makes these objects probably the best SNe 
for this purpose. The relation between the Hubble constant, the 
absolute corrected magnitude M and the zero-point Z in ( [15] ) 
reads 



log 



Ha 



kins 1 Mpcr 



= 0.2 • (M -Z)+5. 



(16) 



Taking the data from SN 1994D, SN 1992 A, and SN 1980N 
presented in Table || and the zero-points from Table || yields a 
Hubble constant of 



H Q = 75±6±6kms- 1 MpcT 1 . 



(17) 



The systematic error originates from the horizontal branch 
magnitude, the universality of the TOM, the adopted Fornax 
distance and the error in Z. The statistical error is just the 
scatter from the three corrected SNe absolute magnitudes. This 
Hubble constant value is dominated by the data from the For- 
nax SNe because of the small statistical weight of SN 1994D. 



5. Summary 

We determined the GCLF distance of NGC 4526 to be 12.0 ± 
3.9 Mpc, placing this galaxy in front of the Virgo cluster. Al- 
though this distance is not precise enough to give a good ab- 
solute magnitude it is sufficient to show that SN 1994D is not 
an overluminous type la SN but rather a normally bright event 
compared to the two Fornax SNe SN 1992A and SN 1980N. 
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All GCLF distances depend on the absolute luminosities 
of the Galactic GCs. It is therefore vital to calibrate the ab- 
solute TOM properly to minimise systematic errors. The most 
important calibration parameter is the absolute magnitude of 
the horizontal branch which directly enters our distance scale. 
The treatment of the Galactic GC sample shows that the asym- 
metric distribution in V and incompleteness effects in the / 
band can introduce a bias in the TOM of over 0.1 mag. 

With the distances and the apparent decline-rates and 
colour-corrected magnitudes of the three SNe SN 1994D, 
SN 1982 A, and SN 1980N one obtains an absolute corrected 
peak magnitude of type la SNe. The Calan/Tololo SN sample 
gives us the zero-point of the Hubble diagram, the second in- 
gredient in calculating the Hubble constant from type la SNe. 
Our data give Hq = 75 ± 6(sys.) ± 6(stat.) kms -1 Mpc -1 . 
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